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ABSTRACT 

Haumea, a rapidly rotating elongated dwarf planet (~ 1500 km in diameter), has two satellites 
and is associated with a "family" of several smaller Kuiper Belt objects (KBOs) in similar orbits. 
All members of the Haumea system share a water ice spectral feature that is distinct from all other 
KBOs. The relative velocities between the Haumea family members arc too small to have formed by 
catastrophic disruption of a large precursor body, which is the process that formed families around 
much smaller asteroids in the Main Belt. Here we show that all of the unusual characteristics of the 
Haumea system are explained by a novel type of giant collision: a graze-and-merge impact between 
two comparably sized bodies. The grazing encounter imparted the high angular momentum that 
spun off fragments from the icy crust of the elongated merged body. The fragments became satellites 
and family members. Giant collision outcomes are extremely sensitive to the impact parameters. 
Compared to the Main Belt, the largest bodies in the Kuiper Belt are more massive and experience 
slower velocity collisions; hence, outcomes of giant collisions are dramatically different between the 
inner and outer solar system. The dwarf planets in the Kuiper Belt record an unexpectedly large 
number of giant collisions, requiring a special dynamical event at the end of solar system formation. 
Subject headings: Kuiper belt objects — individual (Haumea), methods: numerical 



1. INTRODUCTION 



The four largest dwarf planets in the Kuiper Belt form 
a distinct populatio n of bodies with high albedos and 
volatile-rich surfaces ( jSchaller k. Brown 2007 Stansberry 
et al.||2008|). A significant history of collisions is sug- 



gested by the abundance of satellites in this group, which 
is muc h higher than exp ected for the Kuiper Belt as a 
whole ( Brown et al.|2006 ) . Three of the four have known 
satellites: Pluto has three, Haumea (formerly 2003 EL 6 i) 
has two, Eris (2003 UB313) has one, and Makemak e (2005 
FY 9 ) has no substantial satellite ( |Brown| [2008| . The 
size and orbits of these satellites are different from those 
found around smaller (100-km size) Kuiper Belt Objects 
(KBOs). To date, most known satellites around smaller 
KBOs are thought to have formed v ia a still-debated 
capture mechanism (Noll et al. 20081. Hence, a differ- 



ent satellite formation process is needed for the dwarf 
planets, and the most promising mechanism is collisions. 
Recently, numerical simulations support a giant collision 
origin for Pluto's massive satellite, Charon ( Canup 2005 1 . 
However, the formation of the smaller satellites on the 
other dwarf planets has not been studied in detail. 

Haumea, a ~ 1500 km diameter classical belt object 
with a semi-major axis of 43 AU, is a particularly puz- 
zling case as it is also associated with several smaller 
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KBOs with diameters between 70 and 365 km. The 
smaller KBOs share similar orbits and surface properties. 
The associated KBOs have been likened to collisionally- 
produced dynamically and compositionally ass ociated 
"families" t hat are observed in the asteroid belt (Brown 
et al.||2007|). We collectively refer to Haumea, its satel- 



lites and proposed family members as the Haumea sys- 
tem. Haumea has the only known family in the Kuiper 
Belt. The Haumea family members share a deep water 
spectral feature and neutral color dBrown et al.| 2007 
Ragozzine & Brown 2007 Schaller & Brown 2UU8D. The 



water feature is unique m the Kuiper Belt (l-Brown et al. 



n2 
(Br 



20071) and indicative of un usually carbon-free water ice 
( FTmlla-Alonso et al.||2009[) . 

Haumea also has the distinction of being the only 
known highly elongated dwarf planet though its precise 
shape is not known. It h as a spin period of only 3.9 hours 
( Rabinowitz et al.|2006 ), the fastest of all the major and 
dwarf planets. The surface of Haumea is nearly homoge- 
nous wi th the exception of a red spot or fai nt red hemi- 
sphere (Lacerda et al. 2008; Laccrda 2009]); hence, the 
light curve is primarily a reflection of the non-spherical 
shap e. Using the observed li ght curve and rotation pe- 



riod, Rabinowitz et al. (2006) fit a density of 2.6 g cm 



assuming an equilibrium fluid body (a Jacobi ellipsoid). 
Although th e dimensions of Haumea ar e not yet uniquely 
constrained ( Lockwood fc Brown|2009 ), the observations 
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requi re a tri-axial shape (see Table [2]) ( Rabinowitz et al. 
2006 1 . The derived density is greater than the ave rage" 
of ~ 2 g cm" 3 for the largest KBOs ( [Brown|[2008| , al- 
though the density may be smaller with some internal 
friction ( |Holsapple||2007 |. 

Based on the relatively clean water ice surface and 
higher than average bulk density, Brown et al. (2007) ar- 
gue that Haumea is differentiated with a rocky core and 
icy mantle. They propose that the family members and 
satellites are collisionally-derived fragments that origi- 
nated primarily from the icy mantle. The satellites and 
family members are orders of magnitude less massive 
than Haumea, and the family members have a minimum 
velocity at infinity (V qo) of about 150 m s _1 (Ragozzine 



results of the analytic prediction ([2.2 1; and 3) simulate 



& Brown 2007 2009). The velocity dispersion is much. 



less than expected it the Haumea family formed as the 
result of a catastrophic im pact, as in the formation of 

Nesvorny et al.|2006[[MicheT 



asteroid bel t families (e.g. 



et al. 2004). In a catastrophic disruption event, a par 



ent body is disrupted and dispersed such that the largest 
remnant is less than or equal to half the original mass. 
In the gravity regime, the fragments have initial veloc- 
ities relative to the largest remnant comparable t o the 
escape speed (K S c) of the disrupted parent body (Benz 
fc Asphaug||1999 ). Thus, the observed magnitude of the 
velocity dispersion (Voo) of asteroid belt families is sig- 
nificant with respect to V^c of the largest remnant; in 
other words, V esc of the parent body was much greater 
than V csc of the largest remnant in gravity dominated 
disruption events. Unlike most asteroid belt families, 
the velocity dispersion among Haumea family members 
is a small fraction of the escape velocity from Haumea 
(V C sc ~ 900 m s -1 ). Based on the current models of fam- 
ily formation via catastrophic disruption, the Haumea 
family could not have formed by catastrophic disruption 
of a much larger parent body. 
Two impact scenarios have been proposed for the for 



mation of the Haumea system. Brown et al. (20071 
proposed an impact event that falls in tne catstropfnc 
disruption category, which does not agree with the ob- 
served velocity dispersion among fa mily members. To 
explai n the small velocity dispersion, |Schlichting fc San 
( 2009 ) suggest the breakup of a single large moon in or- 
bit around Haumea. However, they do not provide an 
explanation for the initial state: a large moon in close 
orbit around a fast-spinning, elongated planet. To date, 
no known impact scenario explains all of the unusual 
characteristics of the Haumea system. 

In this work, we quantitatively model the formation 
of the Haumea system. We propose that the Haumea 
family formed via a novel type of giant collision: a graze 
and merge impact between two comparably sized bodies 
resulting in high angular momentum, which spun off icy 
fragments that became satellites and family members. 
The analytic and numerical methods are described in |2| 
The results are presented in 33J and the implications for 
giant impacts in the Kuiper Belt are discussed in S|4] 

2. METHODS 

To reduce the parameter space of possible collisions 
that produce a Haumea-like system, we used a three 
step process: 1) derive a n an alytic prediction of plausi- 
ble impact parameters (j |2.1[ ); 2) conduct low- resolution 
simulations over a broad parameter space based on the 



t he m ost promising impact scenarios m high resolution 
( p.3[ ). The next three sections outline the method used 
in each step. 

2.1. Analytic Prediction of Impact Parameters 



Following Canup et al. (2001), using the conservation 
of energy and momentum, we derived an expression for 
the impact parameter and projectile-to-target mass ra- 
tio needed to obtain the observed angular momentum of 
Haumea via a giant impact. Because Haumea is rotating 
near its spin instability limit and the observed velocity 
dispersion among family members is small, we consider 
the case where all of the angular momentum from the 
collision is retained in the remaining body. Assuming 
that the relative velocity between the projectile and tar- 
get was zero at infinity, the impact parameter, b, is given 
by 

b ^ ^ ( f ) 

Lcrit 1/2/(7) 

where b is in units of the sum of the radii of the pro- 
jectile and target, L is the angular momentum, and 
Lcrit = fcM 5 / 3 G 1/2 ( i f^) 1/6 , where k is the inertial con- 
stant (2/5 for a sphere), M is the total mass, and p is 
the bulk density. L cr it is, therefore, the critical spin an- 
gular momentum that a spherical body with constant 
density can sustain. The mass ratio of projectile to 
the total mass, 7 = Mt+m p > en t ers Eq. 1 through 

f{i) = 7(1 - 7)(7 1/3 + (1 - 7 1/3 )) 1/2 sin 6 where is the 
impact angle. However, if we assume that the impact 
velocity is greater than U csc , the to tal energy equation 
will no longer equal zero (Eq. B3 of Canup et al.||2001 ). 
In this case, the impact parameter is constrained by 



b = 



L 



V" 



Lcrit V2/(7) Vi 



(2) 



is the mutual 



where Vi is the impact velocity and V ei 
escape velocity (Eq. 1 of | Canup |[2005| . 

For Haumea we assu me a mass of 4.2 x 10 21 k g and a 
spin period 3.92 hr from |Rabinowitz et al. (2006). Using 
these values, y^— ranges between 1.1 and 0.8 for a plau- 
sible range of bulk densities of the colliding bodies (1.5 
to 2.5 gcm~ 3 ). The results from Eq. Fy and are shown 
in Fig. 

If the impact velocity equals the mutual escape veloc- 
ity of the projectile and target, Eq. [I] requires an impact 
parameter close to one (> 0.8) and a projectile close 
to the mass of the target in order to attain the angular 
momentum in Haumea (Fig. fl]) . Increasing the bulk den- 
sity of the bodies broadens the range of impact parame- 
ters that could transfer the observed angular momentum. 
Similarly, when Vi > V CS c, there is a larger range of pos- 
sible projectile-to-target mass ratios that could produce 
the angular momentum of Haumea (Eq. [2j dotted line in 
Fig. [T] assumes the mean bulk density for large KBOs of 
2 g cm -3 ). This work assumes no initial spin in the tar- 
get and projectile. Initial spin in the same sense as the 
spin angular momentum would potentially increase the 
range of mass ratio that could produce a Haumea-like 
remnant. 
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Fig. 1. — Analytic prediction of impact scenarios leading to the 
angular momentum of Haumea for various bulk densities and Vj. 
Solutions are lines as a function of impact parameter, b, and mass 
of the projectile normalized by the total mass, Mp /Mxot- Eor Vi = 
Vesc (Eq. [lb the short dashed, solid, and long dashed lines show 
solutions for bulk densities of 1.5, 2.0, and 2.5 g cm -3 , respectively. 
For Vi = 1.25V eS c (Eq. [2| the dotted line shows solutions for a 
bulk density of 2.0 g cm " 3 . Haumea parameters are taken from 
|Rabinowitz et al. | p006l l . 

As the impact velocity increases, the remaining body 
does not retain most of the angular momentum of the 
encounter, as is assumed in the equations above. When 
the impact velocity is large enough to begin to disrup t 



the target (3 - 3.5 km s , Stewart & Leinhardt 



2009), 



be car- 



a significant amount of angular momentum will 
ried away by the s mallest fragments (see Fig. 2 in Lein- 



hardt et al.||2000 [). Compensating for the partial loss 
of angular momentum by further increasing the impact 
velocity will lead to the catastrophic disruption regime. 
Recall that the catastrophic disruption regime is ruled 
out based on the velocity dispersion among the family 
members. Although the angular momentum distribution 
in the catastrophic regime has not been extensively stud- 
ied, we expect that angular momentum transfer to the 
largest remnant is inefficient. The largest remnant is a 
gravitationally reaccumulated body; hence, the angular 
momentum of the reaccumulated mass does not approach 
the spin instability limit. Numerical simulations of catas- 
trophic disruption that resolve the shape of the largest 
remnant produce sp herical remnants not fast spinning 
elongated rem nants (Leinhardt et al. 2000 Leinhardt & 
Stewart|[2rjQ9| . 



1'hus, the analytic argument suggests that the most 
straightforward way to create a body with high angu- 
lar momentum via a collision is in a grazing impact be- 
tween two objects of similar size at a velocity close to the 
mutual escape speed. Note, however, that the analytic 
solution makes several assumptions. For example, mass 
that is lost from the system as a result of the collision is 
not taken into account. The escaping mass carries away 
some energy and momentum. In addition, the energy 
and momentum conservation are approximated and do 



not include terms for energy lost in heat and fracturing 
of the target and projectile. 

2.2. Low Resolution Simulations 

To examine the predictions of the idealized analytic 
solutions, we conducted a series of low-resolution simu- 
lations of impacts between gravitational aggregates (Ta- 
ble [I]). Because the encounter velocities are modest 
(< 1.2 km s _1 ), the energy lost to shock deformation 
during the collision is minimal, and gravitational forces 
dominate. Thus, for efficiency, the low-resolution simu- 
lations utilized the TV-body gravity code pkd grav, which 
resolves inelastic particle-particle collisions (|Richardson 
et~aT1[2000{ [Leinhardt et al.|[2000l|Leinhardt & Kichard- 
son||2002l ~^ Particle collisions were modeled using a hard 
sphere model, where the unbreakable spherical particles 
are non-penetrating. The outcome of each inelastic col- 
lision were governed by conservation of momentum and 
typical coefficients of restitution of 0.5 and 1.0 in the nor- 
mal and tangential directions, respectively, for particles 



representing i ce or rock (see discussion in Leinhardt & 
Stewart||2l)09l ) 



Each projectile and target was modeled as a rubble 
pile, a gravitationally b ound aggregate of 955 particles 
with no tensile streng th (Leinhardt et al. 2000[ [Leinhardt 
& Richardson 2002J). Previous simulations ( |Leinhardt 
et al. ' 2000 Fshow that a thousand particles is enough 
to resolve general shape features in a rubble pile. We 
assumed two internal configurations: homogeneous and 
differentiated bodies. In the differentiated cases, the col- 
liding bodies had two layers: a 1.0 g cm" 3 mantle rep- 
resenting ice and a 3.0 g cm -3 core representing rock. 
The mass ratio of the icy mantle to rocky core was also 
varied to reach the desired bulk density, which resulted 
in a range of initial radii from about 500 to 800 km. 

Based on the analytic solutions, the projectile-to- 
target mass ratio was assumed to be one in most cases 
with a subset of simulations considering Mp/Mp =0.5. 
The parameter space included initial bulk densities from 
1.0 to 3.0 g cm" 3 , impact velocities between 0.7 to 1.2 
km 



impact parameters from 0.55 to 0.71, and total 
system masses of 4.5 to 8.2 x 10 21 kg. 

Due to the low resolution, each individual particle had 
relatively large mass, and it was difficult to strip material 
from the surface of the largest remnant. As a result, we 
use these simulations only to refine the impact param- 
eters that reproduce the observed rotation period and 
approximate mass of Haumea. The calculated mass of 
the largest remnant will be biased slightly upward, and 
the family members will be completely unresolved. 

2.3. High Resolution Simulations 

Based on the low-resolution simulation results, impact 
parameters were chosen for high-resolution calculations 
of the formation of Haumea and its family members (Ta- 
ble [2]). We used a hybrid hydrocode to TV-body code 
technique similar to the method used in studies of catas 
trophic collision s in the asteroid be l t and the outer so 
lar system (eg. IMiche l et al.| [20031 JDurda et al 1 12004 



2006; |Lemhardt fe Stewart||2^U^ TrIe 



Nesvorny et al 

hybrid technique captures the shock deformation during 
the early stage of the collision and follows the gravity- 
controlled evolution of the material to very late times. 
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We used GADGET ( |Springel |2005| ) , a smoothed parti- 
cle hydrodynamic s code (SPH) modified to use tabular 



equations of state ( |Marcus et al. 2009), for the hydrocode 
phase and pkdgrav for the gravity pha 



rase of the calcula- 



tion. 



SPH is a Lagrangian technique for solving the hydro- 
dynamic equations in which the mass distribution is rep- 
resented by spherically sym metric overlapping particles 



that are ev olved with time ( Gingold k, Monaghan] 1977 
Lucy|1977 ). SPH has been used extensively to model im 



pacts in the solar system from asteroid collisions and fam 
ily forming e vents to the formation of the Pluto-Charon 
system 



Canup 2005 ) 



Asphaug et aT]|199S 
Although GADGK. 



1998 Michel et al. 2003 



includes self-gravity 



it is not practical to use a hydrocode for the entire inte 
gration of the collision as the timestep is limited by the 
Courant condition. 

The targets and impactors were differentiated bodies 
composed of an ice mantle over a rock core with a bulk 
density of ~ 2 g cm~ 3 . The largest KBOs have bulk den- 
sities of about 2 g cm~ 3 , which is similar to the density 
predicted from cosmochemical es timates of the rock to 
ice ra tio in the outer solar system ( |McKinnon et al.|1997 
2008 1 . Althou gh the internal structu res ot large KBUs 



are unknown (Leinhardt et al. 20081), strong water fea- 



tures on the surfaces ot the largest bodies sugg est that 



they have differentiated (McKinnon et al. 20081. Hence, 



in this work, we consider collisions between differentiated 
bodies only in the high resolution simulations. 

The material in the rocky cores were modeled using a 
tabulated versi on of the mol ecular ANEOS equation of 



state for Si0 2 ( Melosh 
modeled with the tabu 



2007), and the ice mantles were 
ar 5- Phase equation of state for 



H 2 ( |Senft fc Stewart||2008) . The internal temperature 



profile is dependent on the ice to rock ratio and the vis- 
cosity of ice; models indicate that temperatures for a 
Char on size body are gener ally low (e.g., < 300 K after 



4 Ga McKinnon et al.|2008 |. As a result, a constant ini- 
tial temperature of 150 K was chosen. The bodies were 
initialized with hydrostatic pressure profiles. They were 
then allowed to settle in isolation for many dynamical 
times at which point all particles have negligible veloci- 
ties (cms s _1 ) at the specified temperature. The number 
of particles ranged from — 1.2 x 10 5 to — 4 x 10 5 (~ 24-36 
particles per target radius) a sufficient resolution to re- 
solve shock heating and the formation of family members, 
and the results were checked for sensitivity to resolution. 
As in all previous studies of giant im pacts , the materials 
are hydrodynamic (see discussion in £3.3). 

The GADGET simulations were run as long as was 
practically feasible, normally about 60 simulation hours. 
At this time, the collision and mass loss process was 
complete; however, more time was needed to determine 
the long-term orbital stability of material bound to the 
largest object. 

2.3.1. Orbital Evolution of Collision Fragments 

The N-body code pkdgrav was used to integrate the 
orbiting fragments for thousands of spin periods of the 
largest remnant. The GADGET output was translated 
and handed off to pkdgrav. In previous work on as- 
teroid family formation, pairs of pkdgrav particles were 
merged into a single particle after each particle-particle 
collision to reduce computation time (resulting in artifi- 



cial perfectly spherical collision remnants) . In this work, 
because of the significant elongation in the largest rem- 
nant, the shape and rotation rate needed to be preserved 
for the orbital evolution calculation, and particle merg- 
ing could not be used. Hen ce, t he pkdgrav calculation 
utilized inelastic collisions (\ 
and gravitational potential o: 



2.2) to preserve the shape 
the largest remnant. 
However, the number of particles in the largest rem- 
nant in the GADGET simulation (~ 10 5 particles) is 
too large to integrate in pkdgrav because of the compu- 
tational expense of calculating the collisions within the 
largest remnant. Therefore, the largest remnant was de- 
resolved to contain ~ 10 3 particles by placing a grid over 
the body and placing all particles within a grid cell into 
a single particle. Each merged particle had a mass equal 
to the combined mass, a position equal to the center of 
mass position, and velocity vector equal to the center of 
mass velocity. The spin and shape of the largest remnant 
was preserved. The mass of the largest remnant using 
inelastic collisions was compared with a perfect merg- 
ing collision outcome to test the stability of the handoff. 
The mass of the largest remnant was very similar in both 
cases. 

The shape and ice-to-rock ratio of individual smaller 
remnants (the satellites and family members) were not 
resolved in the SPH simulation. Hence, the self- 
gravitating remnants outside of the largest remnant were 
merged into single particles. In this mannner, the de- 
resolved pkdgrav calculation allowed the stability sim- 
ulation to run for thousands of orbits in a reasonable 
amount of time. 

3. RESULTS 

3.1. A Collision Scenario for the Haumea System 

The impact conditions and outcomes of the low- 
resolution numerical simulations of possible Haumea- 
forming impact events are summarized in Table [l] The 
last column of the table, collision type, describes the gen- 
eral class of the collision. We found three different colli- 
sion outcomes in our restricted parameter space (see ^4] 
for further discussion): 1) merge - the projectile and tar- 
get merge after initial impact with little or no mass loss; 
2) graze & merge - the projectile initially hits the target 
with a large impact parameter and then separates, the 
projectile is decelerated, but remains relatively intact, 
and subsequently recollides at a much slower velocity re- 
sulting in a merger and a fast-spinning body; 3) graze & 
run - the projectile and target hit but do not lose enough 
energy to remain bound to each other. 

The impact parameters that produced elongated bod- 
ies with a total mass and spin period similar to Haumea 
tend to be of the graze and merge category. Thus, colli- 
sions that form a Haumea-like body are found in a dis- 
tinctly different parameter space than catastrophic dis- 
ruption events. 

High-resolution hybrid simulations of the most success- 
ful collision scenarios for forming a Haumea-like planet 
were conducted to investigate the properties of the satel- 
lites and family members (Table v3 . A time series from 
an example collision simulation (sim. 4) is presented in 
Fig. [2j which displays the materials in cross-section look- 
ing down on the collision plane. The last frame, which 
has been rotated by ninety degrees to show the collision 
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TABLE 1 

Summary of parameters and results from selected pkdgrav simulations. 



R 


M P 


Pb 




b 


V 


Mi, 


semi-axes of lr 


Pir 




Collision Type 


km 


M T 


g/cm 3 


Pm 


_ 


km/s 


M Tot 


kmxkmxkm 


g/cm 3 


hr 




770 


0.5 


2.0 


1.0 


0.80 


0.7 


0.98 


1270x750x69Cpl 


1.6 


4.4 


merge 


770 


0.5 


2.0 


1.0 


0.80 


0.8 


0.88 


1000 x 708 X67CT 


2.0 


5.0 


graze & merge 


770 


0.5 


2.0 


1.0 


0.80 


0.9 










graze & run 


650 


1.0 


2.0 


1.0 


0.65 


0.8 


0.99 


1396x711x632 


1.7 


4.9 


merge 


650 


1.0 


2.0 


1.0 


0.65 


0.9 


0.99 


1398x750x670 


1.5 


5.0 


graze & merge 


650 


1.0 


2.0 


1.0 


0.65 


1.0 










graze & run 


800 


1.0 


1.0 


1.0 


0.6 


1.0 










graze & run 


650 


1.0 


2.0 


3.0 


0.6 


1.0 










graze & run 


630 


1.0 


2.6 


3.0 


0.6 


1.0 


0.99 


1270x657x607 


2.1 


3.8 


graze & merge 


580 


1.0 


3.3 


3.0 


0.6 


1.0 


0.99 


1164x606x526 


2.8 


3.7 


graze & merge 


650 


1.0 


2.0 


3.0 


0.55 


0.7 


0.99 


1219x717x699 


1.8 


4.1 


merge 


650 


1.0 


2.0 


3.0 


0.55 


0.8 


0.99 


1293x692x666 


1.8 


4.3 


merge 


650 


1.0 


2.0 


3.0 


0.55 


0.9 


0.99 


1385x680x634 


1.8 


4.2 


merge 


650 


1.0 


2.0 


3.0 


0.55 


1.0 


0.99 


1362x708x647 


1.7 


4.5 


graze & merge 


650 


1.0 


2.0 


3.0 


0.55 


1.1 










graze & run 


650 


1.0 


2.0 


3.0 


0.55 


1.2 










graze & run 


650 


1.0 


2.0 


3.0 


0.60 


0.8 


0.99 


1303x702x688 


1.7 


4.1 


merge 


650 


1.0 


2.0 


3.0 


0.60 


0.9 


0.99 


1322x446x636 


1.7 


3.9 


graze & merge 


650 


1.0 


2.0 


3.0 


0.60 


1.0 










graze & run 


650 


1.0 


2.0 


3.0 


0.65 


0.7 


0.99 


1252x740x672 


1.7 


4.1 


merge 


650 


1.0 


2.0 


3.0 


0.65 


0.8 


0.99 


1273x698x697 


1.7 


4.2 


graze & merge 


650 


1.0 


2.0 


3.0 


0.65 


0.9 


0.99 


1445x699x648 


1.6 


4.9 


graze & graze & merge 



a _R - radius of target; - mass of projectile normalized by mass of target; pi,, p c / p m - bulk density, 

density ratio of core to mantle; for pd pm > 1-0, p c = 3.0 and p m = 1.0 g cm" 3 ; b - impact parameter; V 
— first impact velocity; Mi r /AfTot ~~ mass of largest remnant normalized by total mass; pi r - bulk density of 
lr; Pi r — spin period of lr. In all cases, each body contained 955 particles and Mxot = 4.5 X 10 21 kg. Bulk 
density, calculated by circumscribing all particles within an axisymmetric ellipsoid, is always a minimum 
value. 

lr is significantly non-axisymmetric. 

outcome edge-on, shows the surfaces of the largest rem- 
nant and the debris field. During the collision, the rocky 
cores of the progenitor bodies merge, and the resulting 
primary body spins so quickly that it sheds icy mantle 
material from the ends in many small clumps. Some of 
this material is gravitationally bound and some escapes 
from the primary. In this scenario, the satellites and 
family members do not originate from the initial contact; 
instead, they are spun off after the subsequent merger. 
As a result, the Voo of the family members are small; in 
other words, the ejection velocities of the fragments are 
not much greater than the escape velocity of the merged 
primary. 

The analytic and numerical results show that the op- 
timum parameter space to form a Haumea-like planet 
is in an encounter slightly more energetic than merging 
between comparably sized bodies. Some of the merg- 
ing cases that are close to the transition to graze and 
merge also produce an elongated, fast-spinning largest 
remnant; however, from high-resolution simulations, we 
find that such merging simulations eject less mass and, 
thus, do not tend to form families. Although the long- 
term relaxation of the bodies is not considered here, we 
note that impacts with a mass ratio less than one pro- 
duced initially nonsymmetric largest remnants that are 
not consistent with observations. 

The collision parameters that achieve the best agree- 
ment with observations fall in a narrow parameter space: 
an impact speed of 800-900 m s _1 and an impact param- 
eter between 0.6 and 0.65 for equal mass progenitors with 
bulk densities of 2 g cm" 3 . At higher impact velocities 
or higher impact parameters, the two bodies escape from 



each other after the impact. At lower impact velocities 
or lower impact parameters, the two bodies merge and 
less material is ejected as potential family members. Ini- 
tial spin would increase or decrease the impact parameter 
needed to achieve the same total angular momentu m, as 



shown in idealized cases in Leinhardt et al. ( 2000 ) 



3.2. Properties of the Satellites and Family 

The impact conditions defined above reproduce the 
mass, spin period and elongation of Haumea, as well as 
the mass and velocity distributions of the observed satel 
lites and family members. In Fig. [3K, the observed fam 



ily (t riangles and dotted line frorn|Ragozzine fe Brown 



2007 1 are normalized to the number of resolved fragments 



m sim. 4 after 2000 spin orbits to facilitate comparison 
Unresolved fragments have too few particles (N < 10) at 
handoff to be numerically resolved. Haumea and the two 
known satellites are indicated by filled triangles. The di- 
ameter of Haumea family members was derived from the 
absolute magnitude assuming similar albedo to Haumea 
of ~ 0.7, and the diameter was converted into mass as- 
suming a bulk density of 1 g cm~ 3 for a primarily ice 
composition. In Fig. [3)3, the open triangles below zero 
indicate the minimum relative velocities of the known 
family members with respect to the center of the family 



using the recon structed position for Haumea (Ragozzine 
fc Brown|[2007] 



As shown m Fig. [3j3, the total mass of potential satel- 
lites and family members is small: < 0.07Mi r for sim. 1-4 
from Table [2j Although the modeled mass is larger than 
the known satellites and family members (~ 0.01A/i r ), 
one does not expect that all of these bodies have been 
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Fig. 2. — Time series of a graze and merge event: 650-km diameter bodies colliding at 900 m s _1 with an impact parameter of 0.6 
(sim. 4 in Table |2t. Cross section view through the collision plane which is in the page. Field of view is initially 5000x5000 km, increasing 
to 10000X 10000 Kin at 11.1 hours. The last frame (58.9 hrs) shows the system edge on. Color denotes the provenance of the materials: 
icy mantles (cyan and blue) and rocky cores (light and dark grey). For visual clarity in the last frame, the radius of each dot equals one 
fifth of the smoothing length of the corresponding SPH particle, for all other frames the dots are point particles. Although some material 
is exchanged during the first impact, each body remains largely intact after separation. The rocky cores merge after the second impact, 
forming a differentiated primary. The surface of the merged body has distinct patches of ice that originate from each of the precursor 
bodies. The fragments thrown from the merged body are primarily material from the icy mantles. An animation is available in the online 
version of the Journal. 



The Formation of the Haumea System 



7 




Fig. 3. — Mass and velocity distributions for simulations 1-4 in Table[5] (A) Cumulative number of resolved fragments larger than mass 
M versus fragment mass (short dash, dotted, long dash, and filled circles are results from sim. 1-4 after ~ 100 orbits, the open circles are 
results from sim. 4 after ~ 2000 orbits). Triangles and dotted line are the observed family normalized to sim. 4. (B) Cumulative velocity 
distribution of fragments (short dashed, dotted, long dashed, and thin solid line, sim. 1-4, respectively, after 100 orbits, thick black line 
is from sim. 4 after 2000 orbits). The lower absicca is normalized by the escape speed from the largest remnant (model Haumea). The 
left ordinate is the mass of fragments with relative velocities greater than V normalized by the mass of the largest remnant, and the right 
ordinate is normalized by the total mass in smaller fragments in long running sim. 4. Open triangles below zero are minimum relative 
velocities of family members. 



TABLE 2 

Summary of Haumea family forming simulation results and observations. 



Sim. 


M T 


R T 


A' 


b 


Vimpl Vi mp 2 Ml r pl T 


semi-axes of lr 


Pli /h 2 o 


No. 


10 21 kg 


km 






m s~ 1 m s _1 10 21 kg g cm -3 


km X km X km 


hr 



2.25 


650 


1 


0.6 


800 


260 


2.25 


650 


1 


0.6 


800 


240 


2.25 


650 


1 


0.65 


800 


280 


2.25 


650 


1 


0.6 


900 


260 


4.80 


830 


0.22 


0.71 


3000 





Obs^ 



4.3 


2.2 


960 x 870 x 640 


3.6 


0.86 


4.3 


2.1 


1090 x 820 x 680 


3.4 


0.79 


4.3 


2.2 


1100 x 940 x 640 


3.7 


0.73 


4.2 


2.1 


1100 x 810 x 620 


3.9 


0.80 


4.64 


2.3 


1700 x 1500 x 1500 


28 




4.006 


2.6 


1000 x 750 x 500 


3.9 





a Mx - mass of target body; i?x - radius of target body; fi - mass ratio of projectile to target; b - impact 
parameter; Vi mp i — first impact velocity; Vi mp 2 - second impact velocity in graze and merge event; M\ r - mass of 
largest remnant; p\ r - bulk density of lr; P\ r - spin period of lr; /h 2 — i ce mass fraction of all smaller fragments. 
pl r and /h 2 from end of hydrocode phase (50-60 hrs); all other simulation results from end of gravity phase 
(~ 100 orbits). 

Higher resolution version of sim. 1 (N = 4 X 10 5 ). 

c Proposed impact scenario from |Brown et al]|2007|. . 

1 Observed characteristics of Haumea ( jRagozzme ~d Brown||2009| | Rabinowitz et al. 1|2006||. Note: simulations 



presented here attempted to match the mass of 4.21 ± 0.1 X lO^ 1 kg quoted in 



Rabinowitz et al. 



(2006). 
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observed or that all should survive to the present day. 
Hence, the modeled mass of smaller bodies is consis- 
tent with the observed Haumea system. In addition, the 
masses of individual satellite and family members are 
small, with most < 10~ 3 Mi r , which is also in agreement 
with the observed family (Fig. [3k ) . 

Almost all of the family members have a relative ve- 
locity that is < 0. 5V CSC with respect to the primary, with 
80-90% of the mass having a speed < 300 m s _1 and 
35-60% of the mass having a speed less than 150 m s _1 . 
Note that 150 m s _1 is the minimum velocity dispersion 
of the observed family; the true values could be about two 
times lar ger when accounting for the unknown orbital ori- 
entation ( Ragozzine fc Brown|2007 1 . The model predicts 
that negligible mass reaches the typical velocity disper- 
sion among 100-km size KBOs of about 1000 m s _1 . 

As a resolution test one simulation was completed at 
higher resolution (no. 2, dotted line in Fig. [3j N = 4x 10 5 
compared to N = 1.2 x 10 5 for simulations 1, 3 - 5) with 
the same collision parameters as sim. 1. The number and 
velocity dispersion of the remnants are consistent with 
the lower resolution simulations thus we arc confident 
we have reached resolution convergence. 

The modeled satellite and family members are com- 
prised almost entirely of icy mantle material (73-86% by 
mass). The bulk density of the largest remnant increases 
by ~ 10% over the initial density of the progenitor bod- 
ies as a result of the preferential stripping of the lower- 
density ice mantle (Fig. [2]). The modeled bulk density of 
the largest remnant is within the range of uncertainty for 
the density of Haumea. These results explain the shared 
water ice spectral feature of the bodies in the Haumea 
system and the complete lack of bodies with a similar 
water ice spectrum in the general (not dynamically as- 
sociated) KBO population around Haumea. Negligble 
mass was dispersed into the background KBO popula- 
tion in the graze and merge family- forming event. 

In one simulation (no. 4), the orbits of the satellites 
and family members were integrated for 2000 spin orbits 
to assess the stability of the newly formed system around 
the elongated primary. At this time, there are ~ 35 ob- 
jects gravitationally bound and in orbit about the pri- 
mary (the mass in orbit is about 0.0TMi r and the mass 
in family members is about 0.06A/i r ). These bodies have 
eccentricities below 0.9 and instantaneous orbits that do 
not intersect the primary. The bound objects and instan- 
taneous orbital parameters were determined us ing the 



[2005] . Here, 
i evolution of 



companion code (Leinhardt & Richardson 
we do not assess the longer term dynamical 
the collisional system. However, we demonstrate that 
some bodies remain in stable orbits around the elon- 
gated, fast-spinning primar y well after the collis ion event. 



Unlike previously thought (Brown et al. 2006), the for- 
mation of multiple small satellites around Haumea does 
not require accretion in a massive disk. 

3.3. The Effect of Material Strength 

Simulations of giant i mpact events usually utilize a 
SPH code (e.g., this wor k, 
|004 " 



Benz ct al. 1986, 1988; Canup 



2007 



Mannova et al.||20U8[ ). A 

las the advantages of an 



2005| |Benz et al 
agrangian SPH calculation 
arbitrarily large spatial domain and efficiency in track- 
ing small fragments. In an Eulerian (grid-based) code, 
the spatial domain must be decided upon in advance 



and tracking small fragments through the mesh is com- 
putationally very expensive. In previous studies, giant 
impact calculations have neglected material strength on 
the grounds that self-gravitational forces and shock pres- 
sure gradients dominate the problem. Most giant impact 
studies have focused on the end stages of planet for- 
mation and considered hypervelocity impact events that 
generate strong shock waves. Hence, it has been reason- 
able to neglect material strength. For collisions between 
dwarf planets at subsonic velocities, however, it is not 
obvious that strength can be neglected. 

We conducted a few comparison three-dimensional 
simulations using the E ulerian shock physics code CTH 



( McGlaun et al. 1990) with adaptive mesh refinement 
( Crawford 1995]). CTH has the option to include self- 



gravitational f orces using the parallel tree method of 



\ Barnes fc Hut 1986 ). The simulations had a resolution 



ot 31 to 42 km (30 to 40 cells across each initial body). 
In the nominal simulations, each CTH cell is comparable 
in physical size to a single GADGET particle within the 
initial bodies, although the effective resolution in CTH 
is slightly higher because of the differences in smoothing 
lengths between the two codes. Tests at twice the res- 
olution (in each dimension) yielded similar results. The 
CTH calculations utilize the same tabulated equation of 
state models as used in the GADGET SPH simulations. 
Each body was initialized in gravitational equilibrium at 
a constant temperature of 150 K. The radii were 650 km, 
with an ice mantle over a rock core such that the bulk 
density was 2 g cm" 3 . To allow for reasonable calcula- 
tion times, material in cells with a bulk density less than 
0.01 g cm~ 3 was discarded from the calculation. Hence, 
fragments that become the moons and family members 
of Haumea are not modeled in CTH. Only the formation 
of the primary (Haumea) is considered. 

Some simulations were hydrodynamic (no shear 
strength) for direct comparison to GADGET. Other sim- 
ulations utilized a simple friction law (the geological 
yield model in CTH) that represents fr iction in fractured 
(damaged) ice ( Senft k, Stewart||2008[ ) : Y + fiP, where 
Y"o = 0.1 MPa is the cohesion, \x = 0.55 is the friction 
coefficient, and P is pressure. The shear strength is lim- 
ited to a maximum of 0.1 GPa. Shear strength is ther- 
mally degraded as the temperature approaches 273 K. 
The tensile strength was 1.7 MPa. For simplicity, the 
same strength model was used in both the ice and rock 
components 



Collins et al. 


2008) 


crystalline roc 


ks. 



Because the CTH calculations require significantly 
more computational time than the GADGET calcula- 
tions, comparisons between the two codes were made at 
early stages in the impact event. Hydrodynamic CTH 
calculations are similar to the GADGET results (Fig. El. 
The minor differences at late times are due to small dif- 
ferences in the initial conditions: the initial separation 
of the bodies and the better resolution of the ice-rock 
interface in CTH. 

In cases with strength, a number of differences are ob- 
served. The radial oscillations in each body after the ini- 
tial contact are not observed with strength. These oscil- 
lations arise from momentum transfer between the rock 
core and ice mantle; the difference in particle velocities 
from the momentum transfer is damped with strength. 
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Fig. 4. — Time series of sim. 1 in Table [2] using GADGET and CTH with and without strength. View of the cross section through the 
impact plane. Each panel is 4000 X 4000 km. Columns 1-3: Colors represent material (blues - ice mantles; browns - rock cores). Column 
4: Color represents temperature. 
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The number of 'strings of fragments' between the two 
bodies after the initial contact is decreased with strength 
(2.8 hours in Fig. E]). Because of dissipation of energy 
in mechanical worK during the first contact, the second 
contact occurs sooner with strength. The temperature 
increase from the impact is higher with strength because 
of dissipation of mechanical work as heat. Note that 
with and without strength, the cores of the two bodies 
merge quickly after the second contact. The purpose of 
the CTH calculation is not to try to develop a realistic 
strength model for large KBOs, but simply to demon- 
strate that the addition of a reasonable amount of inter- 
nal friction yields essentially the same result as in the 
detailed hydrodynamic simulations. 

3.4. Impact-induced Heating and the Icy Surface 

During the impact event, the temperature in the por- 
tions of the ice mantle in the regions of contact are raised 
to the phase boundary (melting or vaporization, depend- 
ing on the local pressure, Fig. [4|. However, negligible 
mass of ice is melted or vaporized. For the ~ 1 km s _1 
impact velocities considered here and a wide range of 
initial temperatures, ice that experiences the peak shock 
pressu res will be raised to the melting or vaporization 
curve ( |Stewart et aL|2008| ). However, no ice is expected 
to completely melt or vaporize. Thus, the outcome of the 
collision event is not sensitive to the initial temperature 
as neither the rock nor ice components experience sign- 
ficant shock-induced phase changes during the collision 
event. As a confirmation, a test SPH simulation with 
a high temperature rock interior yielded essentially the 
same result. 

In a graze and merge impact, the ice mantles expe- 
rience severe deformation that is expected to fragment 
the ice. Upon merging, solid state differentiation occurs 
quickly, even with some strength. The surface of the re- 
sulting primary (Haumea) is covered with ice fragments 
that have been heated over a wide range: from negligibly 
to the phase boundaries. The fragmented ice mantle has 
negligible cohesion and may be likened to a slurry. Ice 
has exceedingly low friction when warm and at modest 
slip velocit ies which are exceed ed in the scenario consid- 
ered here (Maeno et al. 2003). Thus, we conclude that 
the ice mantle will have negligible resistance to mass loss 
due to the high angular momentum of the merged body. 

Because the graze and merge portions of the impact 
event are similar in the GADGET and CTH calculations 
(with and without strength), we argue that the spin off of 
ice fragments at the end stages of the event is reasonably 
modeled in the hydrodynamic GADGET calculations. 

One of the remarkable characteristics of Haumea, its 
larger satellite, and its family members is the strong wa- 
ter ice spectral feature. The graze and merge family 
forming event explains the ice-dominated surfaces of the 
Haumea system. Although the simulations cannot ad- 
dress the details of ice separation from other phases, the 
results suggest that the formation event produced an icy 
surface t hat is cleaner compared to other KBOs, as is 
observed ( Pinilla-Alonso et al.||200 9). A relatively clean 
icy mantle i s necessary to prevent re ddening from cosmic 
irradiation ( |Rabinowitz et a l. 2008) over the time since 
the impact event (> 1 Ga, Ragozzine & Brown 2007 ) . 
Because the surface of Haumea appears so homogeneous 
and unlike other KBOs, we argue that it is unlikely 



that the surfaces of the precursor bodies were so simi- 
lar. Hence, the strong ice feature on the Haumea family 
must be related to the family-forming impact event. 

3.5. Additional Formation Scenarios 



Brown et al. (20071 suggested a possible impact sce- 
nario to produce the Haumea collisional family based 
on numerical studies of asteroidal family forming events 
( Benz fc Asphaug|T9 99): impact of a projectile with 0.22 



the mass of the target at 3 km s _ at an impact param- 
eter of 0.71. The authors suggested that such an im- 
pact would both strip off a portion of the target's mantle 
and impart a high spin period. We conducted a high- 
resolution simulation of the proposed impact scenario. 
Our results demonstrate that there is not enough energy 
and momentum coupling between the projectile and the 
target to produce a fast-spinning primary (sim. 5 in Ta- 
ble]^. The projectile shaves off some of the target ma- 
terial and escapes from the system, but the remaining 
angular momentum is insufficient to elongate the target 
body. Based on our analytic calculations and numeri- 
cal results, it is not favorable for a small, fast projectile 
to impart enough angular momentum to create a fast- 
spinning elongated body. 

Usi ng an order of magnitude analysis, Schlichting & 



Sari (2009) suggest a different Haumea family forming 



scenario involving a two stage formation process. First, a 
giant impact creates an elongated fast-spinning primary 
and a large, tightly-bound satellite. Second, a subse- 
quent impact onto the satellite, disrupts it and creates 
the family members and small satellites. Based on nu- 
merical simulations, an elongated and fast spinning pri- 
mary is only produced in a slow collision with a large 
impac t parameter and a mas s ratio close to unity (Ta- 
ble [Tj |Leinhardt et al. 2000: L einhardt fc R ichardson 
2002); furthermore, these scenarios do not create a large 
tightly-bound satellite. Impact e vents that pr oduce a 



Canup 20051 do not 



large tightly-bound satellite (e.g. 
form elongated primaries. 

We have not completed an exhaustive parameter space 
study of family- forming collisions in the Kuiper Belt. It 
is possible that additional collisional scenarios such as 
graze and run could form the Haumea collisional family. 
Future work on this scenario would need to address the 
homogenization of the entire surface of Haumea (since 
it would not all melt as a result of the encounter) and 
the probabilities of the encounter providing enough an- 
gular momentum to a Haumea-sized target. In addition, 
a graze and run collision deposits only a small amount 
of internal energy into the target - not enough to dif- 
ferentiate the body. Here we present the first fully self- 
consistent formation model that does not require unusual 
pre-impact conditions. 

4. DISCUSSION 

Giant impacts are common in the late stages of planet 
formation, and several outcomes are possible. With ap- 
proximately decreasing impact energy, a giant impact 
leads to: (i) catastroph ic disruption (half or more of 
the target mass is los t, JBenz fc Asp haug 1999} |Stew-| 
art fc Leinhardt||2009| |Marcus et al.||2009| . The largest 
remnants are aggregates of gravitationally reaccumulated 
material; it is the formation mechanism for asteroid belt 



families, (ii) graze and run (a.k.a. hit and run, Asphaug 



The Formation of the Haumea System 



11 



et al. 2006). In some oblique impacts, the projectile 
hits and then escapes the target with both bodies re- 
maining largely intact. Note that we prefer the term 
"graze" rather than "hit" to indicate the need for an 
oblique impact. (Hi) graze and capture. The projectile 
obliquely hits the target and separates, but does not have 
enough energy to escape and is captured in orbit. It is 
the favored formation mechanism for the Pluto-Charon 

O) 



binary (fCanup 2005|). 
jectile oh 



graze and merge. The pro- 
liquely hits, separates, and then recollides and 
merges with the target. The high angular momentum of 
the merged body spins off some material. Such a sce- 
nario explains the unusual characteristics of the Haumea 
collisional family. In an alternate scenario, a graze and 
merge impact to form Pluto may also p roduce a dis k 
of bound material of the mass of Charon ( |Canup|[2005| . 
(v) merge. The projectile and target merge during the 
first contact; in some cases, a small fraction of material 
may be thrown i nto orbit or escape. The Earth's moon 
forming impact ( Canup|[2004 ) falls in between the graze 
and merge and merging categories. The majority of the 
impact scenarios summarized above have left a distinct 
type of observable satellite or family system in the solar 
system. 

The modeled impact scenarios that reproduce the 
Pluto and Haumea systems are strikingly similar. The 
impact velocities, impact parameters, and masses of the 
projectile are almost exactly the same, with the primary 
difference being the m ass ratio betw een the projectile and 
target: 0.3 for Pluto ( |Canup||2005[ ) and 1.0 for Haumea. 
What is the liklihood of such impact events? The Kuipcr 
Belt is presently composed of multiple dynamically dis 
tinct s ub-populations (e.g., c lassical, scattered and res 



onant, |Morbidelli e t al. 20081. There is strong evidence 
for an excitation event after an initial period of colli- 
sional growth. The Haumea family must have formed 
after the excitation event, which would have scattered 
away all the family members. Although collision proba- 
bilities between two 1000 -km scale bodies in t he classical 
population are negligible, |Levison et al. (2008) calculated 
up to a ~ 50% probability of the impact scenario pro- 
posed by |Brown et al.| ( |2007[ ) if the two bodies originated 
from the scattered population and then entered the clas- 
sical population as a result of the impact. However, the 
estimated mean impact velocity between scattered ob- 
jects is substantially higher (about 2.7 km s _1 ) than we 
find for the Haumea-forming event. The probabilities of 
slower speed collisions need to be investigated. 

5. CONCLUSIONS 

In this work, we show that a graze and merge collision 
event between nearly equal mass bodies is able to pro- 
duce a symmetrically-elongated, fast-spinning primary, 
a family of collisional fragments, and multiple bound 
satellites. The satellites and family members are derived 
from the icy mantle of the merged, differentiated pri- 
mary. The family members have a small relative velocity 
with respect to the primary. This scenario matches all of 
the observed characteristics of the Haumea system. Our 
results predict that observations of future family mem- 
bers should all have relative velocities less than ~ 0.5\4 SC 
of Haumea and that the family members should not be 
isotropically distributed because they formed in a plane. 

We now have the tools to read the record of giant col- 



lisions throughout the solar system. At present, there 
are several IP 's of Pluto-siz ed bodies known in the outer 
solar system ( |Brown|[2008| . The new Pan-STARRS ob- 
servatory will detect up t o an order of magnitude more 
bodies in the Kuiper Belt (Trujillo 2008 1, which will com- 
plete the Haumea system and test the predictions of the 
graze and merge scenario. Models of the dynamical evo- 
lution of the Kuiper Belt indicate that the population 
must have be en much larger in the past (by a factor of 
e.g. ~ 1000) ( jMorbidelli et al.|2008[ ). Satellite formation 
was a common outcome of giant impacts, and the abun- 
dance of satellites around dwarf planets indicates that 
giant impacts were frequent. We find that the types of 
collision events that formed the observed satellites and 
families in the outer solar system are distinctly different 
from the Earth's moon formation and families in the as- 
teroid belt. The narrow range of impact parameters that 
formed the Pluto and Haumea systems place strong con- 
straints on the dynamical history of the largest bodies in 
the Kuiper Belt. 

Acknowledgements. We thank D. Ragozzine and M. 
Cuk for helpful discussions. The hydrocode calculations 
were run on the Odyssey cluster supported by the Har- 
vard FAS Research Computing Group. The N-body cal- 
culations were run using the University of Cambridge, 
Astrophysical Fluids Research Group computational fa- 
cilities. ZML is supported by a STFC fellowship; STS 
by NASA grant # NNX09AP27G. 



12 



Leinhardt, Marcus, and Stewart 



REFERENCES 



Asphaug, E., Agnor, C. B., & Williams, Q. 2006, Nature, 439, 155 
Asphaug, E., Ostro, S. J., Hudson, R. S., Scheeres, D. J., & Benz, 

W. 1998, Nature, 393, 437 
Barnes, J., & Hut, P. 1986, Nature, 324, 446 
Benz, W., Anic, A., Horner, J., & Whitby, J. A. 2007, Space 

Science Reviews, 132, 189 
Benz, W., & Asphaug, E. 1999, Icarus, 142, 5 
Benz, W., Slattery, W. L., & Cameron, A. G. W. 1986, Icarus, 

66, 515 
— . 1988, Icarus, 74, 516 

Brown, M. E. 2008, in The Solar System Beyond Neptune, ed. 

M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, & 

A. Morbidelli, 335-344 
Brown, M. E., Barkume, K. M., Ragozzinc, D., & Schaller, E. L. 

2007, Nature, 446, 294 

Brown, M. E., van Dam, M. A., Bouchez, A. H., Le Mignant, D., 
Campbell, R. D., Chin, J. C. Y., Conrad, A., Hartman, S. K., 
Johansson, E. M., Lafon, R. E., Rabinowitz, D. L., Stomski, 
Jr., P. J., Summers, D. M., Trujillo, C. A., & Wizinowich, P. L. 
2006, ApJL, 639, L43 

Canup, R. M. 2004, Icarus, 168, 433 

— . 2005, Science, 307, 546 

Canup, R. M., Ward, W. R., & Cameron, A. G. W. 2001, Icarus, 
150, 288 

Collins, G. S., Kenkmann, T., Osinski, G. R., & Wunnemann, K. 

2008, Meteoritics and Planetary Science, 43, 1955 
Crawford, D. 1999, Adaptive mesh refinement in CTH, Technical 

report SAND99-1118C, Sandia National Laboratories, 

Albuquerque, NM 
Durda, D. D., Bottke, W. F., Enke, B. L., Merline, W. J., 

Asphaug, E., Richardson, D. C, & Leinhardt, Z. M. 2004, 

Icarus, 170, 243 
Gingold, R. A., & Monaghan, J. J. 1977, MNRAS, 181, 375 
Holsapplc, K. A. 2007, Icarus, 187, 500 
Lacerda, P. 2009, AJ, 137, 3404 

Lacerda, P., Jewitt, D., & Peixinho, N. 2008, AJ, 135, 1749 
Leinhardt, Z. M., & Richardson, D. C. 2002, Icarus, 159, 306 
— . 2005, Icarus, 176, 432 

Leinhardt, Z. M., Richardson, D. C, & Quinn, T. 2000, Icarus, 
146, 133 

Leinhardt, Z. M., & Stewart, S. T. 2009, Icarus, 199, 542 

Leinhardt, Z. M., Stewart, S. T., & Schultz, P. H. 2008, in The 
Solar System Beyond Neptune, ed. M. A. Barucci, 
H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli, 195-211 

Lcvison, H. F., Morbidelli, A., Vokrouhlicky, D., & Bottke, W. F. 
2008, AJ, 136, 1079 

Lockwood, A., & Brown, M. E. 2009, in AAS/Division for 

Planetary Sciences Meeting Abstracts, Vol. 41, AAS/Division 
for Planetary Sciences Meeting Abstracts, 65.06 — h 

Lucy, L. B. 1977, AJ, 82, 1013 

Macno, N., Arakawa, M., Yasutome, A., Mizukami, N., & 
Kanazawa, S. 2003, Canadian Journal of Physics, 81, 241 



Marcus, R. A., Stewart, S. T., Sasselov, D., & Hernquist, L. 2009, 

ApJL, 700, L118 
Marinova, M. M., Aharonson, O., & Asphaug, E. 2008, Nature, 

453, 1216 

McGlaun, J. M., Thompson, S. L., & Elrick, M. G. 1990, Journal 

of Impact Engineering, 10, 351 
McKinnon, W., Prialnik, D., Stern, S., & Coradini, A. 2008, in 

The Solar System Beyond Neptune, ed. M. A. Barucci, 

H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli, 213-241 
McKinnon, W., Simonelli, D., & Schubert, G. 1997, in Pluto and 

Charon, ed. S. Stern & D. Tholen, 295=343 
Melosh, H. J. 2007, Meteoritics and Planetary Science, 42, 2079 
Michel, P., Benz, W., & Richardson, D. C. 2003, Nature, 421, 608 
— . 2004, Planet. Space Sci., 52, 1109 

Morbidelli, A., Levison, H, & Gomes, R. 2008, in The Solar 

System Beyond Neptune, ed. M. A. Barucci, H. Boehnhardt, 

D. P. Cruikshank, & A. Morbidelli, 275-292 
Nesvorny, D., Enke, B. L., Bottke, W. F., Durda, D. D., 

Asphaug, E., & Richardson, D. C. 2006, Icarus, 183, 296 
Noll, K., Grundy, W., Chiang, E., Margot, J., & Kern, S. 2008, in 

The Solar System Beyond Neptune, ed. M. A. Barucci, 

H. Boehnhardt, D. P. Cruikshank, & A. Morbidelli, 345-363 
Pinilla-Alonso, N., Brunetto, R., Licandro, J., Gil-Hutton, R., 

Roush, T. L., & Strazzulla, G. 2009, A&A, 496, 547 
Rabinowitz, D. L., Barkume, K., Brown, M. E., Roe, H., 

Schwartz, M., Tourtellotte, S., & Trujillo, C. 2006, ApJ, 639, 

1238 

Rabinowitz, D. L., Schacfcr, B. E., Schacfcr, M., & Tourtellotte, 

S. W. 2008, AJ, 136, 1502 
Ragozzine, D., & Brown, M. E. 2007, AJ, 134, 2160 
— . 2009, AJ, 137, 4766 

Richardson, D. C, Quinn, T., Stadel, J., & Lake, G. 2000, Icarus, 
143, 45 

Schaller, E. L., & Brown, M. E. 2007, ApJL, 659, L61 
— . 2008, ApJL, 684, L107 

Schlichting, H., & Sari, R. 2009, ApJ, 700, 1242 

Senft, L. E., & Stewart, S. T. 2008, Meteoritics and Planetary 

Science, 43, 1993 
Springel, V. 2005, MNRAS, 364, 1105 

Stansberry, J., Grundy, W., Brown, M., Cruikshank, D., Spencer, 
J., Trilling, D., & Margot, J.-L. 2008, in The Solar System 
Beyond Neptune, ed. M. A. Barucci, H. Boehnhardt, D. P. 
Cruikshank, & A. Morbidelli, 161-179 

Stewart, S., Seiftcr, A., & Obst, A. W. 2008, Geophysical 
Research Letters, 35, L23203 

Stewart, S. T., & Leinhardt, Z. M. 2009, ApJL, 691, L133 

Trujillo, C. A. 2008, in The Solar System Beyond Neptune, ed. 
M. A. Barucci, H. Boehnhardt, D. P. Cruikshank, & 
A. Morbidelli, 573-585 



